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ABSTRACT 

We present a measurement of the R-mode polarization power spectrum (the BB spectrum) from 
100 deg of sky observed with SPTpol, a polarization-sensitive receiver currently installed on the 
South Pole Telescope. The observations used in this work were taken during 2012 and early 2013 
and include data in spectral bands centered at 95 and 150 GHz. We report the BB spectrum in five 
bins in multipole space, spanning the range 300 < i < 2300, and for three spectral combinations: 

95 GHz x 95 GHz, 95 GHz x 150 GHz, and 150 GHz x 150 GHz. We subtract small (< 0.5er in units of 
statistical uncertainty) biases from these spectra and account for the uncertainty in those biases. The 
resulting power spectra are inconsistent with zero power but consistent with predictions for the BB 
spectrum arising from the gravitational lensing of U-mode polarization. If we assume no other source 
of BB power besides lensed B modes, we determine a preference for lensed B modes of 4.9<r. After 
marginalizing over tensor power and foregrounds, namely polarized emission from galactic dust and 
extragalactic sources, this significance is 4.3<r. Fitting for a single parameter, Ai ens , that multiplies 
the predicted lensed R-mode spectrum, and marginalizing over tensor power and foregrounds, we find 
Aiens = 1.08 ±0.26, indicating that our measured spectra are consistent with the signal expected from 
gravitational lensing. The data presented here provide the best measurement to date of the R-mode 
power spectrum on these angular scales. 

Subject headings: cosmic background radiation - cosmology: observations 
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1. INTRODUCTION 

Measurements of the cosmic microwave background 
(CMB), the oldest light in the universe, contain a wealth 
of cosmological information, informing our understand¬ 
ing of physical processes a cross nearly all of cosmic time 
(see Hu & Dodelson 2002 for a review). The majority of 
CMB photons last interacted electromagnetically with 
matter at the epoch of recombination (z ~ 1000), and it 
is that period of cosmic history that is most tightly con¬ 
strained by CMB observations. However, the interaction 
between CMB photons and matter at lower redshifts en¬ 
codes information about the more recent history of the 
universe. In particular, the bending of CMB photon tra¬ 
jectories due to gravitational lensing enables reconstruc¬ 
tion of the gravitational potential between recombination 
and z = 0. Furthermore, the imprint of gravitational 
waves on the polarization of the CMB has the potential 
to probe the absolute earliest moments of cosmic time. 

Current CMB-derived constraints on cosmological pa¬ 
rameters rely primarily on information from the angular 
power spectrum of th e CMB temperatu re fluctuations 
(the TT spectrum, e.g., Planck Collaboration XVI 2014 ). 
However, the polarization of the CMB holds a wealth of 
potential information that is just beginning to be ex¬ 
ploited. As with any headless vector held on the sphere, 
the linear polarization of the CMB 1 can be decomposed 
into a curl-free component and a divergence-free compo¬ 
nent, often called “E modes” and “B modes,” respec¬ 
tively, after the analogous fields in electrodynamics. The 
primary mechanism responsible for CMB polarization is 
Thomson scattering between electrons and CMB plio- 
t ons with an. aniso tropic temperature distribution (e.g., 
Hu & White 1997). Scalar density perturbations at the 
epoch of recombination produce only E-mode polariza¬ 
tion to first order via this mechanism. B-mode polariza¬ 
tion in the CMB can be produced by vector perturbations 
(primordial vorticity) or tensor perturbations (gravita¬ 
tional waves) and by the distortion of E modes through 
gravitational lensing by matt er between the last scatter- 
ing surface and the observer (Seljak & Zaldarriagall 1997 

Seliakpi 


Kamionkowsk i et al.| [l997 Zaidarriaga & Seljak 19118). 


The search for iCmocfe polarization in the CMB has 
been a topic of particular interest because the most suc¬ 
cessful model for explaining many of the observed fea¬ 
tures of the universe, the paradigm of cosmic inflation, 
predicts th e existence of a backgro und of gravitational 
waves (e.g., Abazajian et al. |2015aD . These gravitational 
waves leave their imprint on the CMB through scattering 
at the epoch of recombination (and again at the epoch of 
reionization) through a contribution to the temperature, 
.E-mode polarization, and—most importantly—E-mode 
polarization of the CMB. The gravity-wave contribution 
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1 The circular polarization of th e CMB is observed to be ex¬ 
tremely small, as expected. See e.g. |Mainini et al.[2013| 


to the temperature and .E-mode power spectra is already 
constrained to be too sm all to measure (due to co smic 


variance, see, e.g. Planck Collaboration XVI 2014); 13- 
mode polarization is the only window for measuring this 
signal in the CMB. The amplitude of the gravitational 
wave background is proportional to the expansion rate 
during inflation, and hence to the energy scale of the in¬ 
flationary potential. Thus, a measurement of primordial 
B-mode polarization in the CMB could potentially probe 
physics at energies approaching the Planck scale. 

The CMB B modes induced by gravitational lensing 
of primordial E modes are both an interesting signal in 
their own right and a potential contaminant to the infla¬ 
tionary gravitational-wave (IGW) B modes. In general, 
the signature of lensing in the temperature and polariza¬ 
tion of the CMB can be used to reconstruct the projected 
gravitational potential <f> be tween t he observer an d the 
last scattering surface (e.g., Seljak & Zaidarriaga 1999). 
The reconstructed potential is sensitive to the evolution 
of large-scale structure. In particular, massive neutri¬ 
nos affect the shape of this potential by streaming out 
of small-scale gravitational perturbations and damping 
the growth of structure on these scales. A high-fidelity 
measurement of CMB lensing can in princ iple measure 


the su m of the neutrino masses (e.g., Abazajian et al. 
2015bI. The estimate of <j> from CMB lensing can also be 


used, in concert with a high-signal-to-noise map of the E 
modes, to predict the B-rnode lensing signal. The pre¬ 
dicted B-mode lensing signal can then be cross-co rrelated 
with a direc t B-mode measurement, as in e.g., |Hanson| 
et al.|(|2013) (hereafter H13), or can be used to clean the 
lensing B modes from a direct B-mode measur emen t, 


thereby improving sensitivity t o IGW B modes (|Knox & 
Song||2002| |Kesden et al.||2002[ ). 


Significant experimental progress has been made re¬ 
cently in the Held of CMB B modes. The first detec¬ 
tion of B modes (H13) was m ade in cross-correlation, u s- 
ing CMB data from SPTpol ( Austermann et aT7||2012 )— 
a polarization-sensitive receiver currently i nstalled on 
the 10-mete r South Pole Telescope (SPT, |Carlstrom| 
et al.| 2011)—and cosmic i nfrared background (C1B) 
data from Herschel -SPIRE (Griffin et al. 2010). Simi¬ 
lar cross-correlation measurements have since been made 
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POLARBEAR team also published a measurement of the 


B-mode ang ular power spectrum (POLARBEAR Col- 
laboration|2014a ) in which ~ 2<j evidence tor the l ensing 


signal was seen. Finally, BICEP2 Collaboration (20141 
reported a strong detection of BB power, including a 
component in excess of the ACDM prediction for lensing 
B modes at i ~ 80. The excess has been confirmed in 
yet deeper 150 GHz data on the same area of sky from 


t he Keck A rray experimen t (Keck Array and BICEP2 
Collaborations et al. 2015), but a r ecen t joint an alysis 


dPBlCEP27 Keck, and Planck data (BICEP2/Keck and 
Planck Collaborations et al.|2015 ) has demonstrated that 


at least half of the excess is due to polarized emission 
from galactic dust and that the residual power is consis¬ 
tent with zero IGW signal. Regardless of the interpreta¬ 
tion of the excess B-mode signal, this analysis reported a 
7er detection of lensing B modes at i ~ 200, the tightest 
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direct measurement of lensing B modes to date. 

In this paper, we present a measurement of the BB 
spectrum in the multipole range 300 < £ < 2300, esti¬ 
mated from 100 deg 2 of data taken with SPTpol in 2012 
and 2013. The data used in this work have significant 
overlap with the data used to make the first detection of 
B modes in H13, but there are several key differences in 
the two data sets and analyses. First and most impor¬ 
tantly, the analysis in H13 detected B modes in SPT¬ 
pol data by cross-correlating SPTpol B-mode maps with 
a predicted B-mode template constructed using the E 
modes measured with SPTpol and an estimate of (j> from 
the CIB. Cross-correlation analyses have the attractive 
property that any systematic effect present in only one 
of the data sets will not bias the result. From a sys- 
tematics perspective, the BB spectrum presented here 
is a much more challenging measurement than the H13 
cross-correlation measurement, and the results presented 
here demonstrate that SPTpol cleanly measures the B 
mode power spectrum at angular scales of tens of ar- 
cminutes. Furthermore, the B modes measured in H13 
are necessarily restricted to the lensing signal induced 
by the part of <p traced by the CIB (the CIB is expected 
to have good but not perfect redshift overlap wi th the 
CMB lensing kernel, see, e.g., Holder et al. 2013). The 
BB spectrum presented in this work is sehsitive to all 
contributing signals, including the full lensing signal, the 
IGW signal (though no detection of IGW B modes is ex¬ 
pected in the £ range probed in this work at the current 
levels of sensitivity), and any foreground contamination. 
This work also contains 95 and 150 GHz data from both 
2012 and early 2013, while H13 focused primarily on 150 
GHz data and used only 2012 data. 

This paper is structured as follows: Section 2 describes 
the telescope and receiver. Section 3 describes the obser¬ 
vations used in this analysis. Section 4 details the data 
reduction process from the raw, time-ordered detector 
data to the stage of single-observation maps. Section 5 
describes how we calculate the power spectrum, including 
identifying and correcting for biases in the measured BB 
spectrum, and presents the measured spectrum. We in¬ 
terpret the results in the context of cosmology and other 
B-mode and lensing results in Section 6, and we conclude 
iu Section 7. 

2. TELESCOPE AND RECEIVER 

The SPT is a 10-meter telescope with a wide held of 
view (~ 1 square degree at 150 GHz), designed for con¬ 
ducting large-area surveys of fluctuations in the temper¬ 
ature and polarizat ion of the CMB. The te lescope is de¬ 
scribed in detail in Carlstrom et al. (2011), and f urther 
detail s abou t the optical design can be found in Padin 
etTalTl (f2008). 

After five years of observation with the original SPT- 
SZ receiver, which was sensitive to the CMB intensity but 
not its polarization, the polarization-sensitive SPTpol re¬ 
ceiver was installed in 2012. SPTpol is equipped with 
1536 polarization-sensitive transition edge sensor (TES) 
bolometers, with 1176 detectors at 150 GHz and 360 de¬ 
tectors at 95 GHz. The detectors in the two bands were 
designed an d fabricated independe ntly and are des cribed 
in detail in Henning et al. (2012 150 GHz) and Sayre 
et al. (2012 95GHz). The 150GHz array is composed 
of seven detector modules, each containing 84 pixels, all 


fabricated at the National Institute for Standards and 
Technology Boulder Laboratory. Each module consists 
of a detector array behind a monolithic feedhorn array. 
Incoming power is coupled through the feedhorns to an 
orthomode transducer (OMT), which splits the light into 
two orthogonal polarization states. The 95 GHz array 
consists of 180 individually packaged dual-polarization 
absorber-coupled polarimeters (a total of 360 detectors) 
fabricated at Argonne National Laboratories. Each pair 
of 95 GHz detectors is coupled to the telescope through 
machined contoured feedhorns. The detectors in both 
observing bands are read out using a digital frequency- 
domain multiplexer system with cryogenic superconduct¬ 
ing quantum interference device (SQUID) amplifiers. 

The focal plane is cooled to ~ 250 mK using a commer¬ 
cial pulse-tube cooler and a three-stage helium refrigera¬ 
tor. The TES detectors are then biased to their operating 
point of ~ 500 mK. The entire receiver is maintained at 
~4K, similar to t he SPT-SZ receiver described in |Carl-| 
strom et al. (2011). The secondary optics (including the 
secondary mirror cryostat) are identical to those used for 
SPT-SZ, with the exception of different heat-blocking fil¬ 
ters near the prime focus. Fo r additional details on th e 
SPTpol instrument design, see Austermann et al. (2012). 


3. OBSERVATIONS 

The first eight months of observing with the SPTpol 
receiver (April-October 2012 and April 2013) were dedi¬ 
cated primarily to observations of a 100 deg“ patch of sky 
centered at right ascension 23h30m and declination —55 
degrees. We refer to this field as the SPTpol “lOOd” field 
to distinguish it from the 500 deg 2 survey field, for which 
observations began in May 2013. The 95 and 150 GHz 
data from the first year (2012) of observations of the lOOd 
field were used in H13. The 150 GHz data from 2012 were 
also used to compute the E-mode power spectrum (EE) 
a nd the temperature- E-mode correlat ion spectrum (TE 


Crites et al. (2014 hereafter C14). Story et al. (2014 
used the 150 GHz data from both years of observation of 
this field to reconstruct the CMB lensing potential and 
analyze its power spectrum. The analysis presented here 
is the first to use the 95 and 150 GHz data on this field 
from both years, and the effective white noise levels are 
approximately 17 and 9 /rK-arcmin in polarization at 95 
and 150 GHz. 

Some minor modifications to the instrument were made 
between the 2012 and 2013 observing seasons that could 
affect the analysis in this work. At 150 GHz, one detector 
module was replaced, and the filters used to define the 
band edges were replaced in both bands. Both of these 
modifications most directly affect the shape and overall 
width of the observing bands, which can in turn affect 
beam shape and absolute calibration. For this reason, we 
estimate beams and absolute calibration independently 
for both seasons. 

All observations of the lOOd field used a constant- 
elevation scan strategy, in which the telescope is slewed 
in azimuth back and forth once then stepped a small 
amount in elevation, with the process repeated until the 
full elevation range is covered. The azimuthal scanning 
speed used in all observations was 0.48 degrees per sec¬ 
ond, corresponding to 0.28 degrees per second on the 
sky at the mean elevation of the field, and the constant- 
velocity portion of the scans were between 8.8 and 10.75 
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degrees in azimuth. The elevation step between scans 
was between 13 and 20 arcminutes. 

Only one half of the azimuthal extent of the field is 
observed at one time, in a “lead-trail” strategy that 
allows for ground contamination to be efficiently sub¬ 
tracted w hen ma p s of t he two field halves are differenced 
(e.g., |Pryke et ah [2009). Since no groun d signal is de¬ 
tected in this analysis (see Section |5.8.2|) or in C14, we 
do not use a lead-trail differencing analysis in this work; 
instead, we simply combine each pair of half-field maps 
into a map of the full field. 

We refer to one pass of the telescope, either from left 
to right or from right to left across half of the field, as 
a “scan,” and we refer to a set of scans that cover an 
entire half-field as an “observation.” Each observation 
lasts 30 minutes, and there are a total of roughly 12,000 
observations in the 2012 and 2013 observing seasons, for 
a total of roughly 6000 individual-observation maps of 
the full field. 


size of the pixels used in mapmaking, which is one arc- 
minute for this analysis. A frequency-domain filter with 
very narrow notches at the pulse-tube cooler frequency 
(and the second and third harmonics of this frequency) 
is applied to the data from each detector over the entire 
observation. A total of < 0.01 Hz of bandwidth (< 0.2% 
of the £ x < 6600 band) is removed with this filter, and 
we do not include its effect in the simu latio ns of the filter 
transfer function described in Section 15.11 
The C14 analysis used slightly different filtering 
choices. The scan-direction high-pass filtering was ac¬ 
complished in C14 using polynomial subtraction only (no 
Fourier modes); for this work, we found that the com¬ 
bination of polynomials and Fourier modes resulted in 
a lower level of spurious H-mod e power created in the 
filtering step (see Section 5.5.1 for details). The anti¬ 
aliasing low-pass filter in Cl4 had a cutoff at £ x ~ 10000, 
because the maps used in that analysis had 0.5-arcminute 
pixels. 


4. DATA REDUCTION: TIME-ORDERED DATA TO MAPS 

In their raw form, the time-ordered data for each ob¬ 
servation used in this work consist of one vector of uncal¬ 
ibrated ADC counts for each detector, representing the 
current through the TES as a function of time, and two 
vectors representing the detector pointing as a function 
of time. The first step in processing these data into an¬ 
gular power spectra is to convert the time-ordered data 
into pixelized maps of the Stokes parameters I (or T for 
“temperature”), Q , and U on the sky. The process used 
in this work for making maps from time-ordered data 
closely follows that of C14; we describe the stages of this 
process and point out any differences with C14 below. 

4.1. Time-ordered, Data Filtering 

The data from each detector in each observation are fil¬ 
tered prior to making maps for a number of reasons: to 
remove modes that are strongly affected by low-frequency 
noise, to prevent aliasing of high-frequency noise to lower 
frequencies when the data are binned into map pixels, 
and to eliminate potential contamination from detector 
coupling to the pulse-tube cooler used to cool the re¬ 
ceiver. Data from each detector in each scan are fit to a 
combination of a first-order polynomial (mean and slope) 
and a set of low-order Fourier modes (sines and cosines). 
The best-fit polynomial and Fourier modes are removed, 
resulting in an effective high-pass filter. When maps are 
made from data filtered in this way, the effect of this 
time-domain filtering is a high-pass filter along the scan 
direction (equivalent to right ascension for observations 
from the South Pole) with a cutoff at angular scales of 
roughly three degrees (or an effective multipole number 
in the scan direction of l x ~ 100). 2 To avoid filtering 
artifacts around bright point sources, all sources above 
50 mjy in unpolarized flux at 150 GHz are masked in 
the fitting process. An anti-aliasing low-pass filter is also 
applied to the data from each detector on each scan, 
resulting in a scan-direction low-pass filter in the maps 
with a cutoff of £ x ~ 6600. This filter cutoff is set by the 


4.2. Relative Calibration 

Before the data from all detectors are combined into 
T, Q, and U maps, a factor is applied to the data from 
each detector to equalize the response to astronomical 
signal across each of the two detector arrays (95 and 
150 GHz). The process of measuring and monitoring this 
relative calibration for each array is identical to the pro¬ 
cess described i n C14 an d used in earlier SPT analyses 
(see |Schaffer et ah 2011); we describe the process briefly 
here. 

The relative calibration process involves regular obser¬ 
vations of the galactic HII region RCW38 and of an inter¬ 
nal chopped blackbody source. A 45-minute observation 
of the galactic HII region RCW38 is conducted approxi¬ 
mately once per day (shorter observations for pointing re¬ 
construction are conducted more frequently), while one- 
minute observations of the internal source are conducted 
at least once per hour. An effective temperature for the 
internal source is assigned individually for each detec¬ 
tor (by comparing to the season average of response to 
RCW38). This value times the response of each detector 
to the internal source in the observation nearest a CMB 
field observation is used for relative calibration. We as¬ 
sume that RCW38 is unpolarized, and internal measure¬ 
ments suggest that it is < 1% polarized. In Section 5.9.3 
we discuss the potential for spurious H-mode polariza¬ 
tion caused by small polarization in RCW38. 

Drifts in the internal source temperature are accounted 
for by comparing the average source response across a de¬ 
tector module 3 in each observation to the module average 
over the entire season. Drifts in atmospheric opacity are 
addressed similarly using average RCW38 responses for 
each module. 


4.3. Polarization Calibration 

Before combining the data from all detectors into T, Q , 
and U maps, we must know the polarization properties 
of the detectors. Each detector is designed to be sensitive 
to linearly polarized radiation at a particular angle and 


2 Throughout this work, we use the flat-sky approximation to 
equate multipole number t with 2ir|u|, where u is the Fourier con¬ 
jugate of Cartesian angle on a patch of sky small enough that 
curvature can be neglected. 


3 Although the 95 GHz detectors are fabricated individually, they 
are effectively grouped into four modules of 45 dual-polarization 
pixels (90 detectors) each, based on the wiring of detectors to 
SQUIDs to SQUID readout boards. 
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insensitive to radiation in the orthogonal polarization. 
The two numbers required to characterize each detec¬ 
tor’s polarization performance are the polarization angle 
9 (the angle of linearly polarized radiation at which the 
response of the detector is maximized) and the polar¬ 
ization efficiency r) p (a measure of the ratio of detector 
response to linearly polarized radiation at 9 to radiation 
polarized at 9 + n/2). 

These properties are measured for each detector in ded¬ 
icated observations of a polarized calibration source dur¬ 
ing the Austral summer season. The calibration obser¬ 
vations and the method by which 9 and r] p are derived 
from the data are described in detail in C14. We briefly 
summarize the observations and data reduction here. 

The polarization calibrator consists of a chopped ther¬ 
mal source located behind two wire grid polarizers (one 
at a fixed angle and one that can rotate), physically 
located 3 kilometers from the telescope. One dual¬ 
polarization pixel is pointed at the calibrator, and the 
rotating polarizer is stepped through nearly 180 degrees 
while the response of the two detectors is monitored. We 
fit the response as a function of rotating grid angle to 
a model including 9 and rj p as free parameters. This 
procedure is repeated for all pixels, with multiple mea¬ 
surements per detector where possible. For detectors for 
which the measurements or model fits do not pass data 
quality cuts (~ 25% of the detectors used in this work), 
we assign the median value measured from a subset of 
detectors, namely those that are on the same detector 
module, have the same nominal angle, and did pass data 
quality cuts. We expect this to be a reliable substitution 
since these subsets of detectors were designed to have 
the same angles, and the successfully measured angles 
are consistent with the median angle. 

The median statistical uncertainty on the detector po¬ 
larization angles is 0.5° per detector, and the systematic 
uncertainty on these angles is estimated to be 1°. The 
mean measured polarization efficiency is 98%, and the 
median statistical error on the efficiency is 0.7% per de¬ 
tector. A correlated error in the estimation of all po¬ 
larization angles will result in mixing between E and B 
modes, and this effect is ad dres sed by the cleaning pro¬ 
cedure described in Section l5Al 

4.4. Data Cuts 

We flag and ignore time-ordered data from individual 
detectors on a per-scan basis and a per-observation basis 
using cut criteria that are nearly identical to those used 
in C14. Briefly, we flag data from individual detectors on 
a per-scan basis based on detector noise and the presence 
of discontinuities in the data, including spikes (generally 
attributed to cosmic rays) and sharp changes in DC level 
(generally attributed to changes in the SQUID operating 
point). Data from a particular detector is not used for 
the entire scan if either of these types of discontinuities is 
detected, or if the rms of the data from that detector in 
that scan is greater than 3.5 times the median or less than 
0.25 times the median. The median rms is calculated 
across all detectors on a module. These cuts remove 
roughly five percent of the data. 

Data from individual detectors are flagged at the full- 


observation level based on response to the t wo ty pes 
of calibration observations described in Section 14.21 and 
noise in the frequency band of interest. For a given 
CMB field observation, data from detectors with low 
signal-to-noise response to the closest observation of ei¬ 
ther RCW38 or the internal blackbody source are ex¬ 
cluded from the map for that CMB field observation. 
The noise is calculated for each detector in each obser¬ 
vation by taking the difference between left-going and 
right-going scans, Fourier transforming, and calculating 
the square root of the mean power between 0.3 Hz and 
2 Hz (corresponding to 400 < £ < 2000 in the scan di¬ 
rection, roughly the signal band of interest to this work). 
Data from detectors with abnormally high or low noise 
in this band are cut from an entire observation. We also 
enforce that both detectors in a dual-polarization pixel 
pass all of these cuts; if one detector is flagged, then the 
pixel partner is flagged as well. We find empirically that 
this cut improves the low-frequency noise of the result¬ 
ing maps, presumably by increasing the fidelity of the 
subtraction of unpolarized atmospheric signal. 


4.5. Map-making 

After the cuts described in the previous section are 
applied, the filtered, relatively calibrated, time-ordered 
data from each detector are combined into T, Q 1 and 
U maps using the pointing information, polarization an¬ 
gle, and a weight for each detector. The procedure used 
in this analysis f or polarized map-mak ing i s similar to 
earlier work, e.g., Couchot et al. (19991 and Jones et al. 
(20071. Briefly, for a single observation of the CMB held, 
weights are calculated for each detector based on detector 
polarization efficiency and noise rms—i.e., w oc (rjp/n) 2 , 
where r) p is the polarization efficiency (one value per de¬ 
tector for all observations, based on the p olariz ation cal¬ 
ibration observations described in Section 4.3), and n is 
the noise rms, calculated as described in the previous 
section. The contributions to the weighted estimates of 
the T, Q, and U values for pixel a from detector i are 
then 


= '^ j A tia Wi d ti (1) 

t 

Qioc — ^ ^ Afi a Wj . dti COS 29i 

t 

Uj% = ^2 A tia Wi dti sin 29 i 

t 

where t runs over time samples, du is the value 
recorded by detector i in sample t, w.j is the weight for 
detector i (defined above), A tia is the pointing operator 
encoding when detector i was pointed at pixel a , and 9i 
is the detector polarization angle, and the hat implies 
measured quantities. 

A 3-by-3 matrix representing the T, Q , and U weights 
and the correlations between the three measurements is 
created for each map pixel using this same information. 
The contribution to the weight matrix for the estimate 
of T, Q , and U in pixel a from detector i is 
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WiN, 


IOL 


WiN ia cos 20i 


WiNia cos 2Oi 
WiN ia COS 2 20 i 


WiN^ sin 20i 


WiNi. 


( 2 ) 


where Ni a is the number of samples in which detec¬ 
tor i was pointed at pixel a. The contributions to the 
weighted T, Q , and U estimates and the weight ma¬ 
trix from each detector in a given observing band are 
summed, the weight matrix is inverted for each pixel, and 
the resulting 3-by-3 matrix is applied to the weighted T, 
Q, and U estimates to produce the unweighted map: 

{f,Q,U} a = W- 1 {fw,Q™,uW}. (3) 

We make maps using the oblique Lambert azimuthal 
equal-area projection with 1-arcminute pixels. This sky 
projection introduces small angle distortions, which we 
account for by rotating the Q and U Stokes components 
across the map to maintain a consistent angular coordi¬ 
nate system in this projection. The maps are in units 
of /tKcmb, the equivalent fluctuations of a 2.73 K black- 
body that would produce the measured deviations in in¬ 
tensity. 


4.6. Combining Single-observation Maps 

The cross-spectrunr analysis described in Section |5.3| 
assumes uniform noise properties across all maps used in 
the analysis and all pixels in the maps. Maps made from 
a one half-hour observation (each) of the lead and trail 
halves of the SPTpol lOOd field have sufficiently non- 
uniform coverage that the cross-spectrum analysis per¬ 
formed on these single maps would be significantly sub- 
optimal. We choose to combine the single-observation 
maps into 41 subsets or “bundles” and do th e cross¬ 
spect rum analysis on these bundles. As shown in Polenta 
et al. (2005), the excess variance on a cross-spectrum 
of data equally split N ways as compared to an auto¬ 
spectrum of the full data scales as 1/N 3 , so the excess 
variance with 41 bundles is negligible. 

In a process simil ar t o the detector-level data cuts 
described in Section |4.4[ we calculate the distribution 
of a number of variables over the entire set of single¬ 
observation maps, and we only include in bundles the 
maps that lie within a certain distance of the median 
for all variables. As in C14, we cut maps based on 
the rms noise, the total map weight, the median map 
pixel weight, and the product of median weight times 
the square of the rms noise. In this analysis, we addi¬ 
tionally cut maps based on the number of contributing 
bolometers and on the rms deviation of the time-ordered 
data averaged across all detectors from a model of that 
average time-ordered data (assuming the primary signal 
is from the atmosphere). We include only those maps 
that pass cuts at both 95 and 150 GHz. A total of 4897 
of the roughly 6000 individual observations of the full 
field are included, with approximately 120 observations 
included in each bundle. 


4.7. Absolute Calibration 

The relative calibration process described in Section 
4.2 results in single-observation or bundled T maps that 


are related to the true temperature on the sky by the fil¬ 
ter and beam transfer function and a single overall num¬ 
ber, the absolute temperature calibration. (The Q and 
U maps have an additional factor of the o verall polariza¬ 
tion efficiency, as discussed in Section 4.8 ) The processes 
for estimating the beam and filter transfer function are 
described in Sect ion |5T) we obtain the absolute tempera¬ 
ture calibration using the same method used in C14. We 
briefly describe the method here; for details, see C14. 

The absolute calibration for the 150 GHz SPTpol data 
is derived by comparing full-season coadded temperature 
maps to temperature maps of the same field made previ¬ 
ously with the SPT-SZ receiver, which were in tu rn cal- 

(2013) 

ck 


ibrated through a comparison of the Story et al. 
SPT- SZ temperature p ower spectrum with the 


PTancl 


Collaboration XVI (2014) temperature power spectrum. 
(The uncertainty on the SPT-SZ-to -Planck calibration 
is estimated to be 1.2% in temperature.) The 95 GHz 
SPTpol data are then calibrated by comparing to the 
150 GHz SPTpol data. 

We compare the 150 GHz SPTpol maps used in this 
work to the 150 GHz SPT-SZ maps of the same field by 
creating cross-power spectra, using the same method we 
use for cr oss-s pectra of SPTpol bundle maps, described 
in Section |5.3[ Specifically, we calculate the ratio of the 
cross-spectrum between a full-depth SPTpol map and a 
half-depth SPT-SZ map and the cross-spectrum of two 
half-depth SPTpol maps, corrected for the difference in 
beams between the two experiments. (The SPT-SZ maps 
are made with the same scan strategy and filtering as the 
SPTpol maps, so the filter transfer function divides out 
of the ratio.) We repeat this for many semi-independent 
sets of half-depth SPT-SZ maps and use the distribu¬ 
tion of the cross-spectrunr ratio to calculate the best-fit 
150 GHz absolute calibration and the uncertainty on this 
calibration. A similar process is used to compare 95 GHz 
SPTpol maps to 150 GHz SPTpol maps and to obtain 
the 95 GHz absolute calibration. In the power spectrum 
pipeline, we multiply all SPTpol maps by these calibra¬ 
tion factors before calculating cross-spectra. 

We calculate the absolute temperature calibration sep¬ 
arately for 2012 and 2013 data because the shape and 
overall width of the observing band changed slightly be¬ 
tween the two seasons (see Section [3] for details). In 
all calculations, we treat the 150 GHz SPTpol full- and 
half-depth maps as noise-free compared to 95 GHz SPT¬ 
pol and 150 GHz SPT-SZ full- and half-depth maps, be¬ 
cause the roughly 3x higher noise in the latter two data 
sets dominates the cross-spectrum uncertainty budget. 
The fractional statistical uncertainties on the calibration 
of SPTpol 150 GHz data to SPT-SZ 150 GHz data are 
0.004 for both seasons (2012 and 2013). In quadrature 
with the SPT-SZ-to -Planck uncertainty, this results in 
a total calibration uncertainty of 1.3% at 150 GHz for 
both seasons, highly correlated between the two seasons. 
The statistical uncertainty on the comparison of SPTpol 
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95 GHz data to SPTpol 150 GHz data is < 0.002 for both 
seasons, which contributes negligibly to the total uncer¬ 
tainty budget; thus, the total calibration uncertainty at 
95 GHz is also 1.3% and nearly 100% correlated with 
the 150 GHz uncertainties. All of these uncertainties are 
much smaller than the fractional uncertainty with which 
we expect to measure the BB spectrum, and they will 
make a negligible contribution to the uncertainty on the 
power spectrum or any derived parameters. 


4.8. Absolute Polarization Calibration 

Although the polarization angles and efficiencies are 
measured with a grou nd-based polarization source as de¬ 
scribed in Section |4~3l we allow for additional freedom in 
the normalization of the Q/U maps. This is to account 
for any mechanism that could alter the effective polar¬ 
ization calibration. One such mechanism, the impact 
of electrical crosstalk on our effective polarized beam, 
is discussed in Section 5.9.4 The polarization calibra¬ 
tion is calculated by comparing the EE spectrum mea¬ 
sured in this work to the EE spectrum from the best-fit 
ACDM model for the Pla nck+LENSING+WP+HIGhL 
constr aint in Table 5 of IPlanck Collaboration XVII 
(2014|). The calculation of the measured EE spectrum 
mirrors that of the BB spectrum as described in Sec¬ 
tion J5| The polarization calibration factor for spectral 
bandro € {95 GHz, 150 GHz} is calculated as 


( rym \2 _ \ ^ / ^EE,mxm, theory /^yEE,mxm, data\ / a\ 

V-* cal/ / j f^b J VV 

b 


where b is the Abin index and Wb are simulation-based, 
inverse-variance weights that sum to one. We per¬ 
form this calculation after applying the absolute tem¬ 
perature calibration described in Section |4.7[ We find 
P c 9 a 5 ! = 1.015 ± 0.024 and P ^ 0 = 1.049 ± 0.014, and we 
multipl y the data polarization maps by these factors. In 
Section |5.9| we discuss the impact of the uncertainty of 
these calibration factors on the BB power spectrum mea¬ 
surement. 


4.9. Apodization Mask 

We construct an apodization mask, W, to downweight 
the high-noise regions at the boundary of the SPT¬ 
pol coverage area and to mitigate mode-coupling. The 
apodization mask is constructed as follows. Each map 
bundle has an associated weight array which is approxi¬ 
mately inversely proportional to the square of the white 
noise level at each pixel of the temperature map. We 
smooth each weight array with a a = 10' Gaussian ker¬ 
nel, threshold the smoothed array on 0.05 of its median 
value, and take the intersection of all thresholded ar¬ 
rays across all bundles. We then pad the boundary of 
the intersection map with 5' of zeros and smooth with 
a A = 60' cosine kernel. We perform this process at 95 
and 150 GHz, and the final apodization mask W is the 
product of the individual 95 and 150 GHz masks. 

4.10. Map-space Cleaning 

We perform two map-space cleaning operations on the 
bundle {T, Q,U} maps. These operations are performed 
on both the dat a ma ps and the simulation maps de¬ 
scribed in Section [5Tl 


First, to reduce our sensitivity to bright, emissive 
sources, we interpolate over regions of the map near all 
sources with unpolarized fluxes S 150 > 50 mJy. For each 
source, we replace the values of the map within r < 6' 
of the source with the median value of the map in an 
annulus defined by 6 ' < r < 10 '. 

Second, we filter the maps to reduce our sensitivity to 
scan-synchronous signals. We see evidence for a polar¬ 
ized scan-synchronous signal with ~5 /iK rms at 95 GHz, 
and one with ~0.6 /iK rms at 150 GHz. We clean these 
signals by removing from each map a template that is 
only a function of RA. Due to SPT’s polar location and 
constant-elevation scan strategy, RA is essentially equiv¬ 
alent to the scan coordinate, and any signal that depends 
only on the scan coordinate will depend only on RA. We 
construct the scan-synchronous template for each map 
by averaging the map in bins of RA and then smooth¬ 
ing this binned function with a l°-in-RA Hann function. 
The resulting one-dimensional template is expanded to a 
two-dimensional template and then subtracted from the 
original, two-dimensional map. 

5. POWER SPECTRUM ANALYSIS 

In this section we describe the process by which the 
maps described in the previous section are reduced to 
BB power spectra. 


5.1. Simulations 

A number of steps in the power spectrum analysis rely 
on mock maps, and here we describe the process by which 
the mock maps are generated. 

First, we generate noise-free sky maps that form 
the input to time-ordered-data simulations. The in¬ 
put sky maps are Gaussian realizations of {T,Q,U} 
generate d in the HEALPiX pixelized sphere format 
(iGorski et al.| [2005) with 0.43' (Aside = 8192) resolution. 
The sky maps are generated using {XT, TE, EE, BB} 
CMB power spectra computed u sing the CAMB Boltz¬ 
mann code (Lewis et al. [2000). The ACDM cosmo¬ 
logical parameters input to (JAMB are taken from the 
Planc k+lensing+WP+high L b est-fit model in Ta¬ 
ble 5 of Planck Collaboration XVI (2014). There is no 
tensor power (r = 0), and no foreground power. There 
are two sets of input sky maps: 


• TEB - These sky maps use all of the non-zero 
ACDM CMB power spectra: {TT, TE, EE, BB}. 
The BB spectrum is due to gravitational lensing 
of E-mode polarization. 

• TE - These sky maps are identical to those de¬ 
scribed above, including identical random seeds, 
but do not include E-mode polarization power: 

cp B = 0 . 


We generate 100 realizations of each set of input sky 
maps and convolve them with an azimuthally-symmetric 
beam function. The beam function is measured using 
observations of planets and PKS 2355-534, the brightest 
source in the lOOd field at these observing frequencies. 
We use four distinct beam functions, one for each com¬ 
bination of spectral band and year of observations, to 
generate input sky maps for each spectral band and year 
of observation. 

















Next, we use the SPTpol pointing information to gen¬ 
erate time-ordered data representing mock observations 
of the input sky maps. Before reducing the time-ordered 
data to maps, we simulate the effect of “crosstalk”, the 
mixing of time-ordered data between different detectors. 
We have used observations of the galactic HII region 
RCW38 to measure percent-level crosstalk between some 
detectors. The crosstalk is believed to originate in the 
readout electronics. We define a crosstalk matrix V a b to 
denote the coupling between detectors a and b , and we 
use V ab to mix the simulated time-ordered data of differ¬ 
ent detectors. The main effect of crosstalk for the BB 
analysis presented her e is T —► B leakage, which we dis¬ 
cuss further in Section 5.5.1 After introducing the effect 
of crosstalk, the simulated time-ordered data are reduced 
to bundle maps in a manner that matches the reduction 
of the real data as described in Section [4] 

Finally, we add realistic noise to the simulated bundle 
maps using differences of the data maps. We generate 
a single “realization” of noise for a specific simulated 
map bundle by combining a random half of its ^120 
constituent data maps into one coadded map, the re¬ 
maining half into another coadded map, and differencing 
the two. The resulting difference has no true sky signal 
and should be statistically representative of the noise in 
this particular map bundle, modulo one small difference. 
When making the normal map bundles we combine the 
two halves with their respective weights, but when mak¬ 
ing the noise bundles we combine with equal weights. 
We estimate that this results in the power in the noise 
bundles being ~0.1% larger than the true noise power. 

We generate 200 realizations of noise for each bundle in 
this manner. The constituent data maps are differenced 
differently in each noise realization. In total there are 
200 realizations of noisy map bundles. Each of the 100 
input sky realizations is used twice, while each of the 200 
noise realizations is used once. 


5.2. Constructing E and B 

Here we describe the process by which Ei and Bi, the 
harmonic-space representations of the E-mode and 13- 
mode polarizations, are constructed from the real-space 
{Q, U} maps. First we address small angle distortions 
introduced by our use of the Lambert azimuthal equal- 
area projection. We account for these distortions by ro¬ 
tating the Q and U Stokes components across the map 
to maintain a consistent angular coordinate system. 

To construct the harmonic-space representation of the 
E-mode pola rizat ion we use the standard convention 
( Zaldarriaga||2001 ) 


Ei = Qi cos 2 4>i + U t sin 2<pi 


(5) 


where <f>i is the azimuthal angle of £, Pi = P{PW} is the 
Fourier transform of each apodized map P £ {Q,U }, and 
W is the apodization window. 

To construct the harmonic-space representation of the 
B -mode polarizati on we use the Xb method described in 
Smith & Zaldarriaga (2007). In this method, an interme- 
diate xb map is constructed as the sum of convolutions 
of the unapodized Q and U maps. The kernels for these 
convolutions are compact in real space; each pixel in the 
Xb map is a curl-like linear combination of all Q jU pixels 
within the surrounding square of 5 x 5 pixels. Finally, we 


construct Bi, the harmonic-space representation of the 
B-mode polarization, 

Be = (((- me+m+ 2 )n E{ XB m (6) 

where FJx'bW} is the Fourier transform of the apodized 
Xb map, and the ~ £~ 2 factor accounts for the second- 
order spatial derivative present in the curl-like linear 
combination of the Q and U maps. 


5.3. Cross-spectra 

We employ a cross-spectrum pseudo-6^ approach to es¬ 
timate the BB power spectrum, as described below. The 
first step is to define the two-dimensional cross-spectrum 
between two CMB fields (a, (5) £ {T, E, B}, each coming 
from a spectral band (y m , v n ) £ {95 GHz, 150 GHz}: 


C° 




Re{(a7’T/3^} (7) 


where (i, j) are indices of distinct map bundles. These 
quantities are defined on a two-dimensional Fourier grid 
with resolution S£ = 25. The next step is to calculate 
one-dimensional bandpowers. These are the weighted 
sums of the two-dimensional spectra within a set of l- 
bins {b}, 


Cf r = Cf r H^ n . ( 8 ) 

leb 

The bandpowers are defined on five Abins, each 5£ = 400 
wide, running from An in = 300 to £ max = 2300. The two- 
dimensional Fourier weight H^ ln is defined as 

H r e nn = W[ 2n Z^ n (9) 

where Wl nn is a Wiener filter optimized for detecting a 
ACDM lensed BB power spectrum and ZV m is a binary 
mask. The Wiener filter is defined as 


W pn = 


<3 


B m B n 


TEB ^B m B n 


,TE\ 


Var < C. 


{<r B "' TE } 

^ J SIB 


( 10 ) 


where both the mean and variance are taken across the 
set of 200 noisy simulations. The numerator is the mean 
two-dimensional power spectrum of the signal of interest, 
the lensed BB power spectrum, and the denominator is 
the two-dimensional variance in the absence of this sig¬ 
nal. We note that this filter is optimized to reject the 
hypothesis of no lensed BB power, while the filter opti¬ 
mized to constrain the amplitude of the lensed BB power 
would use TEB simulations in the denominator. Given 
that the data presented here are noise-variance-limited, 
the two filters are quite similar. Both the numerator and 
denominator are smoothed by a Gaussian kernel with 
full-width at half-maximum FWHM; = 150 prior to tak¬ 
ing their ratio. 

The binary mask Z ™ is used to mask out modes satis¬ 
fying \£ x \ < 175 or \£ y \ < 150 for any spectrum involving 
95 GHz data. This mask is necessar y to p ass some of the 
jackknife tests described in Section 5.8. 2| 

We note that while H™ n has been defined to optimize 
sensitivity to the lensed BB spectrum, we simply use the 
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same weight for the other CMB field combinations ( EB , 
EE, etc.), as the BB spectrum is the focus of interest 
for this work. In Figure [l] we show the real-space repre¬ 
sentation of the E-modes analyzed in this work, namely 
B m = Re{ J" -1 { y/Hp™B ^ 1 }}. Each map E m has been 
multiplied by the apodization mask W for the purpose 
of visualization. 


result of this work, the constraint on the amplitude of 
lensing E-mode power, is insensitive to the particular 
choice of k max ; the lensing amplitude shifts by ~0.1er 
when trying & max G {2,4, 6}. 

The linear estimate of the amplitude A of a particular 
form of leakage fft e in £-bin b of the data B™ for spectral 
band m is 


5.4. Convolutional Cleaning 

Given that the CMB temperature and E-mode polar¬ 
ization signals are expected to contain significantly more 
power than the E-mode polarization signal, the potential 
for leakage from T or E into B is a central issue for any 
BB analysis. The leakage from a number of instrumen¬ 
tal effects — gain mismatch, differential pointing, dif¬ 
ferential beam ellipticity, incorrect polarization angles, 
crosstalk, etc. — appears (to first order) as a spatially- 
local leakage of the T or E skies into B. In other words, 
the leakage appears as a convolution of T or E into B. 
Here we describe a process, “convolutional cleaning”, by 
which we reduce the impact of convolutional (T, E) —> B 
leakage. The basic idea is to use the CMB T and E 
modes themselves to estimate and project out convolu¬ 
tional leakage for 0(100) convolution kernels per spectral 
band. We note that, after applying this cleaning, we are 
largely insensitive to any potential parity-violating cos¬ 
mological TB or EB signal. 

We consider leakage of T or E into B that is described 
by a real-space convolution kernel, or equivalently, a mul¬ 
tiplication in two-dimensional Fourier space by a complex 
coupling function ft- For example, the leakage from T 
to B for a particular coupling function ft looks like 


B e = A f t T t . (11) 

The strength of the leakage in the data is set by the pa¬ 
rameter A, and we estimate A in each £-bin b for each 
coupling function. In other words, we make the approx¬ 
imation that the amplitude of a particular form of leak¬ 
age is constant over one £-bin. There is still freedom 
for ^-dependence of the leakage because the amplitude is 
allowed to vary between £-bins. 

We use coupling functions of the form 

ft 9 = cos (kept + 9) i k (12) 

where k is a small integer k G {0, 1,..., fc max }, 4>e is 
the azimuthal angle in two-dimensional Fourier space, 
9 G {0, 7t/2}, and the i k factor ensures that / corresponds 
to a purely real (rather than imaginary) real-space ker¬ 
nel. We use this form because the coupling functions 
for well known types of leakage can be expressed as lin¬ 
ear combinations of couplings with the form described 
above. For example, differential gain has a TB coupling 
ft oc cos (2 </><), differential pointing has & TB coupling 
ft oc ii (cos(<^/) + cos( 3</>£)), and a global polarization 
angle offset has an EB coupling ft oc 1. We emphasize 
that we do not assume a particular physical mechanism 
(such as differential pointing or differential ellipticity) 
for this leakage. Rather, we assume that the leakage is 
varying slowly in £, and has an azimuthal dependence 
described by cos (kept + 9) for small integer k. We use 
fc max = 4, which allows for the cleaning of any of the 
types of leakage described above. We find that the main 


A; 


kGXr. 


Zte b Re {(B?rXtf k “}H t 


jmm 


E, eb wm. 


mm 

£ 


(13) 


where X G {T, E}, and the weight Hff 1 ™ 1 is defined in Sec¬ 


tion 5.3 The numerator is linear in the data B7 1 , and the 
denominator is used to normalize the estimate. Because 
the 150 GHz data are deeper than the 95 GHz data, we 
use the 150 GHz Tt and Et modes to estimate the leak¬ 
age in the 95 and 150 GHz Bt modes. We use a cross¬ 
spectrum approach to estimate both Re{(B]p)*Xtf]f 0 } 
and \Xtfi 6 \ 2 . For our baseline value of fc max = 4, there 
are a total of 18 leakage modes per £-bin, or 90 modes 
for all five £-bins, per spectral band. 

Ideally we would use these estimates of the leakage am¬ 
plitudes to remove from the data B™ the inferred leaked 
E-mode, 


E 


leak 


’ m = E 4 


kOXm 
b 


Xtf k t e q b t 


(14) 


kOXb 


where the two-dimensional binning operator q^t is one 
for t G b and zero otherwise. In practice, however, due 
to the non-orthogonality of the different modes being re¬ 
moved, the cleaning process does not converge after one 
iteration. We address this issue using a damped, itera¬ 
tive scheme. In each iteration we estimate the amplitude 
of each form of leakage; subtract a damped version of the 
inferred leaked mode, 0.2E^ cak ’ m , from Bp-, and repeat 
for all forms of leakage. We iterate over these steps 20 
times, at which point the process has converged. We use 
the accumulated estimate for each leakage amplitude to 
construct the final inferred leakage mode, Ei eak,m . Fi¬ 
nally, we subtract B l p ak,m from E^ 1 ’* for each bundle i. 

To be clear, the same Ei eak,m is subtracted from all bun¬ 
dles. In other words, all two-dimensional Bt that appear 
in this work have been cleaned in this way, with one ex¬ 
ception: the Fourier weights Ht are calculated using Bt 
that have not been cleaned. 

We find that the convolutional cleaning process 
strongly suppresses the one form of convolutional leak¬ 
age that we expect to be present in the SPTpol data: 
crosstalk-induced leakage. We have used our simulated 
maps to determine that convolutional cleaning reduces 
the additive bias to the EE spectrum introduced by 
crosstalk by factors of ~20-100. The cost of this im¬ 
provement, however, is the introduction of a small, ad¬ 
ditive noise bias, as discussed below in Section 5.5.1 and 
in the Appendix. 


5.5. Unbiased Spectra 

Here we describe the process by which the raw EE 
bandpowers described in Section |5.3 are processed into 
unbiased bandpowers. The raw bandpowers are subject 
to additive and multiplicative biases, and we correct for 
these biases using simulated bandpowers. 
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Fig. 1. The real-space representation of the 5-modes analyzed in this work, at 9 5 GH z (top) and 150 GHz (bottom). The 5-modes 
have been filtered by the two-dimensional Fourier weights y/Hg, as defined in Section |5.3| and multiplied by the apodization mask. The 
crosshatch pattern faintly visible in the 95 GHz map is due to being zero near £ x = U and £ y =0. These maps are noise-dominated on 
all angular scales. 
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5.5.1. Additive Bias 

Our strategy for dealing with additive bias — any BB 
power that is present in the absence of a true B-mode po¬ 
larization signal — is straightforward. We measure the 
mean BB power present in noisy, simulated bandpowers 
generated using TE input skies (i.e., no true BB power), 
and we subtract this bias from the data bandpowers and 
the simulated TEB bandpowers. This subtraction is per¬ 
formed prior to correcting for any multiplicative bias. 
The additive biases can to some degree be organized as 
follows. 

• E —► B from geometry and filtering - Imperfect 
separation of E and Bona small area of sky and 
the high-pass filtering of the time-ordered data re¬ 
sult in E —> B leakage. The resulting bias in BB 
is approximately +0.6er in the lowest Abin of the 
150 x 150 spectrum, and smaller than 0.1a in other 
Abins of this spectrum or any Abin of the other 
spectra. 

• (T, E) —> B from crosstalk - The electrical 
crosstalk between detectors results in (T, E) —> B 
leakage, predominantly T —> B. Although the raw 
crosstalk-induced additive bias is as high as 1.0a in 
the lowest Abin of 150x150, convolutional cleaning 
reduces this to < 0.1a. 

• Negative noise bias from convolutional cleaning - 
We find that the convolutional cleaning step intro¬ 
duces a negative bias. By analyzing two sets of sim¬ 
ulations which have convolutional cleaning turned 
on and off, we can isolate this effect. We find 
that its magnitude and shape are in broad agree¬ 
ment with a simplified analytic calculation given 
in the Appendix. The amplitude of the bias is ap¬ 
proximately —0.4a in the lowest two Abins of the 
95 x 95 and 150 x 150 spectra, smaller elsewhere, 
and approximately zero at 95 x 150. The convolu¬ 
tional cleaning process has strongly attenuated bi¬ 
ases arising from convolutional leakage (crosstalk, 
beam systematics, etc.) in exchange for a small 
noise bias that we can characterize effectively per¬ 
fectly. 



Fig. 2.— The net additive bias to the three BB spectra as a frac¬ 
tion of the diagonal statistical uncertainty. Both of these quantities 
are determined using noisy simulated data. The biases are < 0.5cr 
and are subtracted from the measured spectra. 

application of the apodization mask (W). K can be ex¬ 
panded as 

K w = P M ( M U '[W ] F e )Q ev . (16) 

where Fp is the effective transfer function from the beam 
and the filtering of the time-ordered data, Php is the bin - 
ning operator and Qev its reciprocal |Hivon ~et al.|2002 ). 
The mode coupling kernel Mpp> [W] accounts for the mix- 
ing of power between bins due to the apodization mask. 
The mode coupling kernel is calculated analytically, as 
described in the appendix of C14. The calculation cor¬ 
rects only for B —»• B coupling and effectively assumes 
that there is no E —> B coupling. The latter is of course 
not true, but we correct for E —► B leakage by subtract¬ 
ing off an additive bias, as described in Section [5.5.1| As 
in C14, the effective transfer function Fp is calculated 
by comparing the mean simulated spectra to the theory 
spectrum input to the simulations. The calculation is 
iterative and converges after two iterations. 


The additive biases are shown in Figure [2] We subtract 
these biases from the measured spectra, and we address 
systematic uncertainties associated with the additive bi¬ 
ases in Section [5791 

5.5.2. Multiplicative Bias 

After removing the additive bias we correct for the mul¬ 
tiplicative bias. The bandpowers Cb are a biased estimate 
of the true binned sky power, CA, due to effects such as 
time-ordered data filtering, beam smoothing, finite sky 
coverage, and mode-mode mixing from the source and 
apodization mask. The biased and unbiased estimates 
are related by 


pB m B n 


= K B m B r 

— n-bb' 


cB™ 3 " 


(15) 


where the K matrix accounts for the effects of the in¬ 
strumental beam and time-ordered data filtering and the 


5.6. Bandpower Covariance 

We approximate the covariance between BB bandpow¬ 
ers as completely diagonal, and we take the variance of 
bandpowers from the set of 200 noisy simulations as our 
estimate of the diagonal variances. There are two sets 
of covariances: those that use TE input skies and thus 
naturally account for the noise variance and the variance 
from T and E leakage, and those that use TEB input 
skies, which additionally account for variance from true 
BB power. The TE variances are useful for rejecting 
the C BB = 0 hypothesis. The TEB variances are only 
slightly larger than the TE variances, at most 20% larger, 
implying that noise power dominates over signal power, 
as expected. 

We have used the noisy simulated bandpowers to as¬ 
sess the validity of the diagonal covariance approxima¬ 
tion. We attempted to measure the covariance between 
off-diagonal spectral combinations e.g., (95 GHz x 95 
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GHz) and (95 GHz x 150 GHz), or between neighboring 
Gbins. In each case, the distribution of covariance esti¬ 
mates is consistent with zero. This is not surprising given 
that the covariance is dominated by noise power and the 
noise is expected to be uncorrelated between spectral 
bands and between £-bins. We have also explicitly tested 
the importance of the off-diagonal covariance terms at 
(95x95)x(95x150) and (150xl50)x(95xl50). We esti¬ 
mated these terms using an analytic approximation, al¬ 
tered the covariance matrix accordingly, and found that 
the main result of this work, the constraint on the ampli¬ 
tude of lensing B-mode power, changed by ~0.1<r. For 
simplicity, we do not include these off-diagonal terms. 

We note that our simulation input skies contain Gaus¬ 
sian realizations of lensed BB power and therefore do not 
account for the ^20% inter-^-bin correlation p resent for 


non-Gaussian, truly-lensed B-modes (Benoit-Levy et al. 


2012|. Given that the non-Gaussian structure leads to a 


20% bin-bin correlation on a source of power (lensed BB 
power) that itself contributes only 10-20% of the band- 
power variance, we ignore the non-Gaussian contribution 
in our TEB variances. 

Finally we note that, because our simulations were free 
of foreground power, we have slightly underestimated the 
BB variance. Power from randomly distributed, polar¬ 
ized point sources should be negligible; even the upper 
limit of the Poisson powers considered in Section [b] are 
< 1% of the data noise power. Variance from polar¬ 
ized galactic dust should also be small compared to the 
noise p ower . For our nominal dust model described in 
Section 6.1 the dust power is approximately 1% of the 
noise power in the lowest £-bin of the 150 GHz x 150 
GHz spectrum, and smaller elsewhere. 


5.7. Bandpower Window Functions 

We compute bandpower window functions to allow 
the measured bandpowers to be compared to theoretical 
power spectra. The window functions, w\/t, are defined 
through the relation 


C b = {w\/ i)C t . 

(17) 

Following the formalism described previously 
tion|5.5.2[ this can be rewritten as 

in Sec- 

C b = (K~ 1 ) bb f P b i a 

(18) 

which implies that 


w\/£ = (K~ 1 ) bb > Pvt' MjP(F(. 

(19) 

5.8. Null Tests 



In this section we test the internal consistency of the 
data by considering two types of null tests. First we 
consider the TB and EB spectra. Second we consider a 
suite of “jackknife” tests that are sensitive to potential 
instrumental systematics. 

5.8.1. TB and EB spectra 

Although we are interested primarily in the BB power 
spectrum, here we consider the other spectra that contain 
the B-mode polarization field, namely the TB and EB 
spectra. We do not expect a cosmological signal in these 
spectra. 


Typically one would compare the TB and EB data 
spectra to the null hypothesis, namely C TB = C EB = 0. 
Here, due to the added complexity of the convolutional 
cleaning step, we do not necessarily expect the simulated 
TB and EB spectra to have zero mean, and indeed we 
find that several of the simulated bandpowers have nega¬ 
tive mean values. Nevertheless, we can test whether the 
data spectra are consistent with the distribution of simu¬ 
lated spectra. We find that the TB and EB data spectra 
at 95 x 95, 95 x 150, and 150 x 95 are all consistent with 
the distribution of simulations. However, the TB and 
EB data spectra at 150 x 150 are systematically lower 
than the mean of the simulated spectra. TB is low by 
[-2.3, -2.1, -1.8, -1.4, -1.4]er in the different Gbins, while 
EB is low by [-1.5, -1.5, -3.0, -0.6, -1.5]cr. 

While we cannot offer an explanation for this difference 
between the data and the simulations, we argue that it 
is not significant for the main focus of this work, the 
BB spectrum. More specifically, we use the difference 
between the data and simulation spectra to estimate the 
resulting spurious BB, and find that it is small compared 
to our statistical uncertainties. For example, for TB the 
difference between the data and simulations is 5C BB = 
^TS,data_^TS,sims^ an( } ^he est i mate 0 f associated 

spurious BB power is C BB = (SC EB ) 2 /C BT . We find 
the spurious BB estimated from either TB or EB to be 
very small, less than 0.2% of the expected ACDM lensed 
BB spectrum. 


5.8.2. Jackknives 

We perform a suite of “jackknife” tests to further val¬ 
idate the consistency of the data. In each jackknife test 
the bundle maps are sorted according to a metric de¬ 
signed to trace a potential systematic effect, and “jack¬ 
knife maps” are constructed by differencing high-metric 
bundle maps with low-metric bundle maps. The resulting 
jackknife maps should not contain sky signal and should 
have a power spectrum, the jackknife spectrum, consis¬ 
tent with noise. We construct the jackknife spectrum 
by taking the mean cross-spectrum among the jackknife 
bundle maps. We note that we do not apply the convo¬ 
lutional cleaning step to the bundle maps going into the 
jackknife maps. The convolutional cleaning step would 
have no effect, as it would remove the same modes iden¬ 
tically from each of the maps prior to taking their differ¬ 
ence. 

We perform four jackknife tests. 

1. Left-Right: Jackknife maps are made by differ¬ 
encing data from left-going and right-going scans. 
This tests for any power or systematic effect that 
is present more strongly in one scan direction, such 
as spurious power caused by the step in telescope 
elevation at the end of each left-right scan pair. 

2. Ground: Jackknife maps are made by dividing the 
maps into two subsets thought to have different 
susceptibility to ground contamination. We use 
the same azimuthal range metric used in SPT-SZ 
power spectrum analyses such as Slrirokoff et al.l 
|2011 This tests for spurious power from nearby 
buildings or other sources of contamination that 
are fixed to the ground. 
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3. 1st lialf-2nd half: Jackknife maps are made by dif¬ 
ferencing data from the first and second halves 
of the set of chronologically ordered map bun¬ 
dles. This tests for any systematic associated with 
the changes made to the SPTpol receiver between 
2012 and 2013, systematics associated with small 
changes to the scan strategy made in late 2012, or 
any other slowly-varying systematic. 

4. Visual Inspection: Jackknife maps are made by dif¬ 
ferencing map bundles that show visually-identified 
anomalies (e.g., faint stripes or spots) and those 
that do not. This tests for systematic power asso¬ 
ciated with these features. 

We use a slightly different notation for the jack¬ 
knife bandpowers than is used for the main bandpow- 
ers. Each jackknife bandpower is , where / £ 
{BB, EB,TB} denotes a CMB field combination, s £ 
{95x95,95x150,150x150} denotes a spectral combina¬ 
tion, and j £ {LR, GROUND, TIME, VIS} denotes a 
jackknife. 

We test the consistency of the jackknife spectra with 
the null hypothesis as follows. We estimate the diagonal 

covariance of the jackknife bandpowers, a 2 us¬ 

ing the variance among the cross-spectra divided by the 
number of unique cross-spectra. Next, we define jack¬ 
knife “x bandpowers”, 



and calculate four test statistics to determine the com¬ 
patibility of this set of spectra with the null hypothesis. 
These test statistics are: 

• max/sj G J^bX-b^ l) * This tests for spectra which 

are preferentially positive or negative across the full 
multipole range. 

• max/sj (j - This tests for spectra which 
preferentially have outlying bandpowers. 

• ma x bfsj ((; xl sj ) 2 ) - This tests for any particularly 
strong outlying bandpower. 

• J2bf SJ (x f b Sj ) 2 - This tests for a general tendency to 
have outlying bandpowers. 

We compare the data values of these test statistics 
to those obtained using a set of 10000 zero-mean, unit- 
width Gaussian realizations of each • We calculate 
the probability to exceed (PTE) the value of each data 
test statistic given the values in the set of simulated 
test statistics. Finally, we define a global test statis¬ 
tic, Pjointj the probability to simultaneously exceed all of 
the test statistics, and calculate the probability to exceed 
1 — Pjointi again using simulations. 

As the focus of this work is the BB spectrum, our nom¬ 
inal set of jackknife CMB field combinations is simply 
{BB}. We find that the PTEs for the four test statis¬ 
tics are 0.29, 0.60, 0.15, and 0.09, and the global PTE is 


0.22. We take this as evidence that the BB spectra do 
not have significant contamination under these jackknife 
tests. 

We have also repeated these jackknife tests using an ex¬ 
panded set of CMB field combinations, {BB,EB,TB}. 
Due to the high signal-to-noise imaging of T and E 
modes, these tests are in principle more difficult to pass. 
In this case, the PTEs for the four test statistics are 0.41, 
0.59, 0.39, 0.02, and the global PTE is 0.14. We take 
this as further evidence that the P-mode data used in 
this work does not have significant contamination under 
these jackknife tests. 

5.9. Systematic Uncertainties 

Here we discuss several potential sources of systematic 
uncertainty in the BB power spectrum measurement. In 
all cases we demonstrate that the systematic uncertainty 
is much smaller than the statistical uncertainties and can 
be safely ignored. 

5.9.1. Uncertainty in Additive BB Bias 

As described in Section |5.5.1[ we use simulations to 
determine the additive biases which must be subtracted 
from the measured BB bandpowers. Here we assess the 
accuracy with which we have determined these biases. 

First we consider the uncertainty due to using a finite 
number of simulation realizations. We have used 200 
realizations, leading to uncertainties on the mean biases 
of 0.07cr per Abin. 

Next we specifically consider T —i B leakage, which is 
dominated by crosstalk leakage. While the raw bias is 
as high as l.Ocr in the lowest Pbin of 150x150, convolu¬ 
tional cleaning reduces this bias to at most O.lcr per t-bin. 
The uncertainty in the absolute temperature calibration 
is less than 3% in power, resulting in a negligible 0.0030- 
uncertainty in the additive bias. Similarly, the underly¬ 
ing ACDM TT spectrum is constrained to approximately 
the same level of precision, and its uncertainty can be 
safely ignored. 

Next we consider E —► B leakage, which is caused pri¬ 
marily by the field geometry and the filtering of time- 
ordered data. This leakage results in a bias as high as 
0 .6(7. The uncertainty in the absolute polarization cali¬ 
bration is less than 5% in power, resulting in a negligible 
0.03(7 uncertainty in the additive bias. Again, the un¬ 
derlying ACDM EE spectrum is constrained to a similar 
level of precision. 

Finally we consider the negative noise bias from con¬ 
volutional cleaning. This bias is approximately -0.4er per 
£-bin in the 95 x 95 or 150 x 150 spectra. However, be¬ 
cause we use the data to generate the noise realizations 
used in the simulations, there is essentially zero system¬ 
atic uncertainty in the level of this bias. 

5.9.2. Uncertainty in Multiplicative Bias in BB 

We discussed above how the uncertainty in the abso¬ 
lute polarization calibration results in a (negligible) un¬ 
certainty in the additive BB bias from E -A B. It also 
results in an uncertainty on the amplitude of the mea¬ 
sured BB spectrum. The uncertainty is 3% and 4% in 
power for the 150x150 and 95x150 spectra respectively, 
where nearly all of our sensitivity lies. This results in 
a ~0.1(7 global systematic uncertainty on the amplitude 
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of the BB spectrum. We note that this small uncer¬ 
tainty does not affect the significance with which we de¬ 
tect BB power, as the statistical uncertainties are noise- 
dominated, and the noise would suffer from the same 
mis-calibration. 

5.9.3. T —> Q/U leakage 

We see evidence for 0.65% ± 0.15% leakage of T into Q 
and U at 150 GHz. We measure this using TQ and TU 
spectra as described in C14. The source of this leakage 
is not known, although it could arise from our celestial 
calibration source, RCW38, being ^0.5% polarized. Un¬ 
like previous SPTpol analyses which explicitly cleaned 
this leakage from the Q and U maps, this work relies on 
convolutional cleaning to remove this leakage. We expect 
the convolutional cleaning process to reduce the associ¬ 
ated additive BB bias by at least a factor of 20, resulting 
in a bias that is at most ~0.15cr in the lowest Gbin of 
150x150 and smaller elsewhere. The leakage is smaller 
yet at 95 GHz. 

5.9.4. Effect of Crosstalk on Beams 

Because the electrical crosstalk between detectors hap¬ 
pens to be preferentially negative, the array-averaged 
temperature beam will tend to lose solid angle. The 
array-averaged polarization (Q/U) beams, on the other 
hand, will not lose solid angle, due to partial cancellation 
of crosstalk from sets of nearly randomly-oriented detec¬ 
tors. This results in the array-averaged temperature and 
polarization beams differing slightly. Throughout this 
work we have used the effective temperature beam as 
measured with observations of planets and AGN. This 
implies that our polarized power spectra have been de- 
biased using slightly incorrect beam functions, and that 
our simulations were performed using slightly incorrect 
beams. The ratio of the effective temperature and po¬ 
larization beam functions can be broken down into an 
G independent mean offset and an ^-dependent shape 
around that mean offset. The ^-independent offset is 
naturally accounted for when we calibrate our TT and 
EE spectra independently. The ^-dependent shape is 
not accounted for, but is well approximated by a linear 
tilt from +2% to -2% in power across the i range of this 
work. This results in negligible (< O.lcr) systematic un¬ 
certainties in the removal of additive and multiplicative 
BB biases. 

5.9.5. Time-dependent Crosstalk 

In our baseline simulations, the crosstalk matrix V a b 
that encodes the coupling between detectors a and b 
was fixed. In fact, we have evidence that there is some 
amount of time-variation in the crosstalk matrix. Time- 
varying crosstalk could potentially introduce a different 
level of additive bias in the BB spectrum than time- 
independent crosstalk does. We have addressed this issue 
with a second set of simulations in which the crosstalk 
matrix V a b was allowed to vary, per observation, in a way 
that mimics the time-variation we observe. We find that 
the resulting additive BB bias is slightly smaller than in 
the time-independent case. The difference is 0.2(7 in the 
lowest G bin of the 150 x 150 spectrum, and smaller than 
0.04u in all other Gbins. 


5.9.6. Small-scale Beam Features 

We have used observations of Venus to measure the 
array-averaged beams on small scales (r < 3'). The re¬ 
sulting polarization maps can be used to place upper 
limits on the T —► B leakage due to anomalous beam 
features on these small scales. We estimate the spurious 

j-) j-) / fyTB ,Venus Venus\9 g _ yT'T'?-A.CDM! /oa ’\ht 

BB power as (C i ’ /C e ) C £ ’ /20. We 

have scaled the ACDM TT spectrum by the appropriate 
beam factor and have divided by 20, the minimum level 
of improvement that convolutional cleaning provides on 
simulated crosstalk-induced leakage (a stand-in for other 
types of convolutional leakage). We find that the result¬ 
ing spurious BB would introduce a net bias less than 
0.02cr. 

While the Venus maps provide an upper limit on the 
T —i B leakage from small-scale beam features, they do 
not address E —>• B leakage. However, given that the 
limits on T —i B are so small, and given that EE is a 
factor of ~50 smaller than TT at these multipoles, we 
do not expect E —>• B from small-scale beam features to 
be significant. 

5.9.7. Large-scale Beam Features 

We use observations of the moon to estimate that ap¬ 
proximately 5% of the beam solid angle is contained in 
the radial range 3 ' < r < 25'. If we consider a pes¬ 
simistic scenario in which this portion of the beam, when 
averaged over all detectors, couples T —► B with 5% effi¬ 
ciency or E —> B with 50% efficiency (both in map units), 
and that convolutional cleaning reduces the resulting BB 
power by a factor of 20, the residual BB power is less 
than 0.03cr for all Gbins. 

5.9.8. Systematics Summary 

To summarize, we have considered a number of sys¬ 
tematic uncertainties and demonstrated that they are 
small compared to the statistical uncertainties. When 
added in quadrature with the statistical uncertainties, 
the systematic uncertainties associated with additive bi¬ 
ases would increase our total uncertainty by ~3% in the 
lowest Gbin of the 150x150 spectrum, and would have 
a smaller effect elsewhere. The systematic uncertainty 
in the global multiplicative bias is ~3.5%. As discussed 
in Section [6j this is significantly smaller than the ^20% 
precision with which we measure any BB signal. Ad¬ 
ditionally, this multiplicative uncertainty does not affect 
the significance with which we detect a BB signal, as 
it would affect the bandpowers and their uncertainties 
nearly equally. 

5.10. Bandpowers 

The final BB bandpowers are provided in Table [l] and 
shown in Figure [3] The bottom panel of Figure |3] shows 
the inverse-variance-weighted combination of the three 
sets of bandpowers. The spectrally-combined bandpow¬ 
ers are for visualization purposes only; the likelihood em¬ 
ployed in Section[6]uses the original set of three bandpow¬ 
ers. The bandpowers, covariance matrix, and bandpower 
window functions are available at the SPT website 4 . 

4 http://pole.uchicago.edu/public/data/keislerl5/ 
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Fig. 3.- Top: BB power spectrum bandpowers from the individual 95GHzx95GHz, 95 GHzx 150 GHz, and 150 GHzX 150 GHz spectra. 
Bottom: The inverse-variance-weighted combination of the three sets of bandpowe rs in the top panel. For reference, the expected lensed 
BB spectrum from the Planck+lensing+WP+highL best-fit model in Table 5 of|Planck Collaboration XVI|(|2014| is shown by a solid 
gray line in each plot. 
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TABLE 1 

BB BANDPOWERS, ^center Cl [10 _3 //K 2 ] 


•^center 

£ range 

95 x 95 
tCi a{lCt) 

95 x 150 
tCi a{iC t ) 

150 x 150 

ICi c(iC t ) 

Combined 

ICt a(£C e ) 

500 

300-700 

3.4 

1.4 

0.88 

0.55 

0.57 

0.33 

0.76 

0.28 

900 

700-1100 

2.9 

1.1 

1.18 

0.50 

0.51 

0.32 

0.82 

0.26 

1300 

1100-1500 

0.4 

1.3 

0.27 

0.50 

1.07 

0.37 

0.77 

0.29 

1700 

1500-1900 

0.2 

1.5 

-0.23 

0.55 

0.16 

0.38 

0.04 

0.30 

2100 

1900-2300 

-1.7 

1.6 

-0.59 

0.61 

-0.29 

0.47 

-0.47 

0.36 


BB bandpowers and uncertainties measured in this work. The last two columns give results for the inverse-variance-weighted combination 
of the three sets of bandpowers. 
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Fig. 4.— BB power spectrum m easurements from SP Tpol (this work), ACTpol (Naess ct ah 2014), BICEP2/Keck (|Keck Array and| 
BICEP2 Collaborations et al.|2015|), and POLARBEAR ([POLARBEAR Collaboration 2U14hj). The Highest multipole bin ot tne Auipol 
data is no t shown. T he solid gray line shows the expected iensed BB spectrum trom the BLANCK+lensing+WP+highL best-fit model in 
Table 5 ofjPlanck Collaboration XVI|(|2014|. The dotted line shows the nominal 150 GHz BB power spectrum of Galactic dust e mission 
used in this work. Tins m odel is derived trom an analysis of polarized dust emission in the BICEP2/Keck field using Planck data (|Planck| 
[Collaboration et al.|20l4) . The dash-dotted line shows the sum of the Iensed BB power and dust BB power. 
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6. INTERPRETATION 

6.1. Significance of BB leasing measurement 

The bandpowers from this work shown in Figure [3] 
show a clear preference for non-zero (positive) power. 
Using a multivariate Gaussian likelihood a nd th e band- 
power covariance matrix defined in Section |5.6[ we find 
that the x 2 of the data to the null hypothesis of zero BB 
power is 41.5, with an associated PTE for 15 degrees of 
freedom (dof) of 3 x 10 -4 . The data are not well fit 
by zero sky power. A far better fit is the predicted lens- 
ing BB spectrum from (for ins tance) t he Planck+ L ens 
best-fit model from Table II of IPlanck Collaboration XVII 
(2014). With zero free parameters, the x 2 of our data to 
this model is 18.1 for 15 dof, for a PTE of 0.25, and a 
Ay 2 relative to zero BB power of 23.4. 

To translate this preference into a detection signifi¬ 
cance, we introduce a single parameter, Ai ens , which we 
use to artificially scale the predicted lensing BB power 
from our fiducial cosmological model. We explore this pa¬ 
rameter space (and all parameter spaces discussed sub¬ 
sequently in this work) using the Markov chain Monte 
Carl o method provided by the CosmoMC 5 software pack¬ 
age (Lewis & Bridle|2002 ). The x 2 for the best-fit model 
is 17.6 (PIE = 0.23). The data prefer the addition of 
this single parameter by a Ax 2 value of 23.9, translat¬ 
ing to a 4.9cr detection of lensing, under the assumption 
of no other sky components. Note that all of these x 2 
values are calculated with noise variance only (i.e., no 
sample variance from sky signals), because we are only 
asking at what level the data in this patch of sky prefer a 
component that looks like lensed B-modes, and are not 
attempting to relate the amplitude of this component to 
any global cosmological parameter. 

We do expect other sky components to contribute to 
the measured BB power, so we also fit for the amplitude 
of the lensing B modes in the presence of other signals. 
We know there will be contributions from polarized emis¬ 
sion from extragalactic sources and from dust emission 
within our own Galaxy, and we add parameters to the fit 
describing each of these. 

The area of sky used in this work is contained within 
the BICEP2 area, and should thus have a similar level 
of galactic polarized du st emission . M otivated by th e 
recent Planck results ( Planck Collaboration et al.|[2014), 
we model the contribution from Galactic dust as: 

Df st (^i x u 2 ) = x v 2 ) = (21) 

Z7T 


Aiust -#80 (^o X Z'o) 


$dust iyi )*Sd ust O2) ( i 


-SdustM 


80 


-0.42 


where C'/ ust (iq x vfi) is the contribution of dust 
to the 95 GHz x 95 GHz, 95GHzxl50GHz, or 
150 GHzx 150 GHz spectrum; Adust is an overall 
scaling of the amplitude of the dust power spectrum 
at all frequencies; = 150GHz; g st (i/o x v$) is the 
best-fit value of the dust BB spectrum in the BICEP2 
observing regio n at 150 GHz and l = 80, according 
to Planck Collabor atio n et al. (2014); and S’dust is 


the Planck Collaboration et al. (2014) assumption for 


temperature). We integrate the dust spectrum over the 
measured SPTpol bandpasses. The only free parameter 
in our fit to this model is Adust- For reference, the 
nominal value of Adust = 1 corresponds to powers of 
Dgust = (0.00169, 0.00447, 0.0118) /iK 2 in the SPTpol 
spectra at (95 x 95, 95 x 150, and 150 x 150). 

We model the contribution from extragalactic sources 
as a Cf = constant term, as would be expected if the 
emission was dominated by the Poisson noise in the num¬ 
ber of sources in our observing region, rather than the an¬ 
gular clustering of those sources. We allow the amplitude 
of the Poisson term in our three frequency combinations 
to vary independently, giving us three free parameters, 
Aps, 95x95, Aps,95x150; an d Aps. 1 50x150- 


Based on t he arguments in BIC EP2 Collabor atio n 
and BICEP2/Keck and Planck Collaborations 


the spectral behavior of the dust - (In units of CMB 


et al. (2015), we expect the contribution from Galactic 
synchroiron emission to be far below our ability to detect 
it, and we ignore this contribution in our analysis. We 
also ignore the contribution from the polarized emission 
by clustered extragalactic sources. We expect this sig¬ 
nal to be far less important relative to the Poisson signal 
in polarization measurements than it is in temperature 
measurements, because the clustered signal is dominated 
by dusty, star-forming galaxies (DSFGs), which we ex¬ 
pect to have a much smaller polarization fraction than 
the synchrotron-emitting active galactic nuclei (AGN) 
expe cted to contribute the bulk of the Po isson signal 
(e.g., Seiffert et al.||2007 Battye et al.||2011 ). 

Finally, we add a component to our lit that has the 
shape of the expected IGW B-mode signal, and we use 
the free parameter r to control the the amplitude of 
this component. We do not expect to be able to distin¬ 
guish this component from Galactic dust in the t range 
and sensitivity level in this work; we include the IGW 
component primarily to facilitate direct comparison be¬ 
tween our l ensing results and those from BICEP2 (see 
Section |6.2| below) . 

Our nominal fit includes priors on all five nuisance pa¬ 
rameters (Aps,95x95; Aps,95xl50; Aps, 150x150; Adust, and 
r). The nominal prior on Adust is a Gaussian centered 
on 1.0, with a = 0.3. The nominal priors on Poisson 
power are uniform between zero and four times the value 
calc ulated for each param eter by int egrating source mod- 
els (jDe Zotti et al. 2010 for AGN, Negrello et al. 2007 
for DSFGs) up to the unpolarized flux cut of 50mJy 
and assuming polarization fractions of 5% for AGN and 
2% fo r DSFGs (the upper bounds from |Seiffert et al. 
2007 and Battye et al.|20lf respectively). For reference, 
the upper edges of Poisson Ct> priors are 1.9 x 10~ 7 ^K 2 , 
9.1 x 10 _8 /zK 2 , and 4.4 x 10" 8 /xK 2 at 95x95, 95x150, 
and 150x150, respectively. The prior on r is uniform be¬ 
tween 0 and 0.4. We also place a prior on Ai ens (uniform 
between 0 and 3.0), but the posterior on Ai ens in all fits 
is dominated by the data, not the prior. 

The best-fit point in this nominal six-parameter space 
(including priors) has a x 2 value of 17.9. If we fix Ai ens = 
0 and re-fit, the best-fit point has a y 2 value of 36.8. 
Thus, the data show a preference of Ay 2 = 18.9, or 4.3 ct, 
for lensing B modes when we marginalize over foreground 
parameters. 


5 http://cosmologist.info/cosmomc/ 
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6 .2. Best-fit lensing amplitude and comparison with 
previous results 

To report a value of A lens that can be used to assess 
the validity of the assumed cosmological model and to 
compare to other values in the literature, we repeat the 
six-parameter fit including BB sample variance in the 
uncertainty budget. The best-fit value and ler uncer¬ 
tainty from this fit is 


A lpns = 1.08 ±0.26 


( 22 ) 


consistent with the value of 1.0 that we would expect 
if the cosmological model we assumed were correct. The 
constraint on Ai ens is not strongly dependent on the fore¬ 
ground priors: if we increase the upper limit of the Pois¬ 
son priors to 100 times the nominal power, the best-fit 
value of Ai ens shifts by less than 1 <j, and the error bar in¬ 
creases by less than 10%. Similar behavior is seen when 
the Adust prior is changed to a uniform prior between 0 
and 3. 

Figure [3] demonstrates that the measurement of BB 


which other sources of 13-modes dominate — is rejected 
at > 3cr independent of whether or not we marginalize 
over foregrounds. 

Next we compare to measurements of the lensing 13- 
mode power spectrum that rely on cross-correlation with 
tracers of the CMB lensing potential </>. In these analy¬ 
ses, a lensing 13-mode template is constructed by lensing 
measured E modes by an estimate of (j> derived from CIB 
maps or CMB lensing. The 13-mode template is then cor¬ 
related with the measured B modes to estimate the lens¬ 
ing BB power spectrum. The linear relationship between 
the CIB and (j> is, in turn, based on the C/~ CIB spectrum 
measured by Planck , and it would therefore be surprising 
if these CIB cross-correlation analyses gave results that 
were not consistent with the (jxj) spectrum measured by 
Planck (which itself is consistent with the ACDM pre¬ 
diction). The best-fit value of Ai ens from the SPTpol 
CIB cross-correlation analysis in H13 was 1.092 ± 0.141. 
The POLARBEAR collaboration does not report a value 
of Aiens from their CIB cross-correlation analysis. The 
ACTPol collaboration measured Ai ens = 1.30 ± 0.40 in 


power in this work is at least visually consistent with 

their cross-correlation analysis with Planck CIB data 

previous measurements in the same i range. The most 


van Engelen et al.| 2014 

)■ 

Planck Collaboration et al. 

straightforward way to compare the results quantita- 


2015C also detect lensed B 

modes in this fashion and 


comparison is complicated somewhat by the fact that 
these values are in general reported with respect to BB 
spectra predicted using different values of cosmological 
parameters. The differences between the predicted am¬ 
plitudes are typically much smaller than the lcr uncer¬ 
tainty on Aiens from any of the measurements, however, 
and we will ignore them in the following comparison. 

First we compare to other i3-mode power spec¬ 
trum measurements. The best-fit value of Ai ens from 
POLARBEAR measurements of the i3-m ode p o wer spec- 


the CIB as a <f> tracer or Ai ens = 0.93 ± 0.08 when using 
the quadratic-estimator-derived <f>. These measurements 
are consistent with the ACDM prediction and with our 
measured value. 

Finally, we note that our constraint on Ai ens is also 
consistent with determinations of this parameter from 
quadratic-estimator reconstructions of the lensing po¬ 
tential using the four-point function of CMB temper¬ 
ature and pol arization data. The strongest such con- 


trum was 1.12±0.61 (statistical only, POLARBEAR Col¬ 
laboration 1 2014a), consistent with the value we find here. 
The ACTPol collaboration does not repor t a val ue o f 
Aicns fro m their 33-mode power s pectr um (Naess et al. 
2014). BICEP2 Collaboration (20141 reports a 5.50- 
detection of lensed 33 modes and a best fit value of 
Ai ens ~ 1.75, roughly 2er above 1.0, with no marginaliza¬ 
tion over foregrounds. We have used the publicly avail¬ 
able BICEP2 likelihood module 6 to repeat this analysis 
with a marginalization over foreground and IGW ten¬ 
sor power that mirrors the SPTpol analysis presented 
here. We obtain a constraint from the BICEP2 data of 
Aiens = 1.45 ± 0. 38, ro ughly 1 .2tr ab ove 1.0. Mo re re¬ 
cently, BICEP2/Keck and Planck Collaborations et al. 
(2015) reported a strong (7er) detection of lensing B- 
rnodes at i ~ 200: Ai ens = 1.13 ± 0.18. This constraint 
was obtained after marginalizing over contributions from 
IGW tensor power and polarized galactic dust emission, 
and represents the most precise direct measurement of 
lensing 33-modes to date. To summarize, these BB power 
spectrum measurements are consistent with each other, 
and with the ACDM prediction. These measurements 
suggest that the BB spectra in these particular fields and 
at these observing frequencies are dominated by lensing 
i3-modes, at least at i > 200. The SPTpol BB spectrum 
presented here is particularly useful in that regard, in the 
sense that Ai ens >2 — a rough proxy for a scenario in 

6 http://bicepkeck.org 


straint (Planck Collaboration et al. 2015) uses full- 
mission Planck data to yield Ai ens = 0.983 ± 0.025. 
Additionally, three recent works use p olarizati on-only 
quadratic estimators to constrain Ai ons . POLARBEAR 


Collabora tion! (|2014c|) find Ai ens = 1.06 ± 0.47 (statisti¬ 


cal onfy)7 Story et al.| ( |2014 ) use SPTpol data to measure 
Aiens = 0.92±0.25 (statistical only ) on the same 100 deg 2 
field u sed in this work, and |Planck Co llaborati on et aL] 
(20151 find Aiens = 1.252±0.350. Again, these measure- 
ments of Ai ens are consistent with the ACDM prediction 
and with our measured value. 

7. CONCLUSION 

We have presented a measurement of the 33-mode 
power spectrum (3333 spectrum) using data from 
100 deg 2 of sky observed with SPTpol in 2012 and early 
2013. The 3333 spectrum is estimated in the multi¬ 
pole range 300 < i < 2300 for three spectral com¬ 
binations: 95 GHz x 95 GHz, 95 GHz x 150 GHz, and 
150 GHz x 150 GHz. These data provide the best mea¬ 
surement of the 33-mode power spectrum on these angu¬ 
lar scales to date. 

Several sources of bias— all at a level below the sta¬ 
tistical uncertainty in the power spectrum—are identi¬ 
fied and subtracted from the data. The resulting power 
spectrum is strongly inconsistent with zero power but 
consistent with predictions for the 3333 spectrum aris¬ 
ing from the gravitational lensing of 33-mode polariza¬ 
tion. In a six-parameter fit that includes the predicted 
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lensed B-mode spectrum scaled by a single parameter 
Ai ens , as well as contributions from Galactic dust, extra- 
galactic sources, and any IGW B-mode signal, we find 
aliens = 1-08 ± 0.26. The null hypothesis of no lensed 
B-modes is ruled out at 4.3cr after marginalizing over 
foreground parameters (4.9cr if foregrounds are fixed to 
zero). 

Improved constraints on the BB spectrum are ex¬ 
pected soon from a number of ongoing CMB experi¬ 
ments. For example, in December 2014 SPTpol com¬ 
pleted the second of three years of observation of a 500 
deg 2 field, and the resulting data will significantly im¬ 
prove upon the BB spectrum presented in this work. 
A future generation of instruments aims to further con¬ 
strain inflationary B modes, provide maps of lensing B 
modes over large fractions of the sky, and constrain or 
measure the mass in neutrinos. We can expect significant 
progress in this field in the coming decade. 
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APPENDIX 

We provide here a simplified calculation of the negative, additive BB bias caused by the convolutional cleaning 
process. We demonstrate that the bias is negative, with a magnitude and shape similar to what we measure using 
simulations. To be clear, we correct for the bias using simulations, not the results of this calculation. The following 
arguments apply for (95 GHz x 95 GHz) and (150 GHz xl50 GHz) spectra. We have no evidence for and do not 
expect a negative noise bias at (95 GHz x 150 GHz), because the noise is uncorrelated between the two spectral bands. 

Consider the uniformly-weighted BB bandpower at Gbin b: 


r*BB 

b 


1 


BbTl pairs 


EE R e{(^r s i} 


( 1 ) 


where i is the two-dimensional Fourier space vector, rib is the number of two-dimensional Fourier grid points that 
belong in Gbin 5, i and j are map bundle indices, and ?r pa irs is the number of map bundle pairs. Note that, for 
purposes of clarity, the notation in the main text does not include the factor of 1 /rib, while we must explicitly include 
that factor here. 

We consider an artificial example in which the B maps contain only noise. In this limit, the bundle pairs are 
uncorrelated and the expectation value of the cross-spectrum is zero: 

(G b BB )cx(^Re{(^)*Hi})=0. (2) 


We will now demonstrate that, after applying the convolutional cleaning process, the expectation of this cross-spectrum 
is negative, with a magnitude similar to the bias we measure using simulations. 

The cleaning process for one convolution kernel ft in one Gbin b essentially amounts to 1) constructing _B leak , the 


projection of the coadded B( onto a unit-length mode gt = 


Xlh 
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and 2) subtracting B ^ eak from each bundle B\ 


The cleaned BB cross-spectrum is then 
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( 3 ) 


( 4 ) 


( 5 ) 


The expectation value of the first term (oc (Bf)*Bj) is still zero. It is straightforward to show that the expectation 
values of the second and third terms are N/nb and —2 N/rib, respectively, where N is the noise power level in the 
coadded B map (assumed here to be Gindependent, i.e., white), and rib is the number of two-dimensional Fourier grid 
points belonging to Gbin b. The cleaned cross-spectrum is then 


(jBB, clean _ 

b n b n b n b 


( 6 ) 


The difference in the expectation values of the cleaned and original power spectra, namely (— N/rib — 0) = —N/rib, is 
the negative additive bias caused by projecting out a single convolution kernel ft. When multiple, orthogonal kernels 
are projected out, the bias will scale as nkemeis- The final bias is then 


(jBB, bias,clean _ jy ^kernels 

b n b 


( 7 ) 
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We compare the results of this calculation to those obtained using the more realistic, simulated bandpowers. We 
assume white noise levels of 17 and 9 /jK-arcmin at 95 and 150 GHz, rikemeis = 18, and rib — /sky (^5 max — C'l min ). 
The resulting biases are somewhat smaller than those obtained using simulations (tlreory/sims ~ 0.6), but the sign, 
shape, and overall magnitude of the two methods are in broad agreement. This demonstrates that the basic mechanism 
of the bias can be understood using this simplified calculation, despite ignoring details such as non-white noise and 
non-uniform Fourier weights. 



